Spatially homogeneous ground state of the two-dimensional Hubbard model 
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We investigate the stability with respect to phase separation or charge density-wave formation of 
the two-dimensional Hubbard model for various values of the local Coulomb repulsion and electron 
densities using Green-function Monte Carlo techniques. The well known sign problem is particularly 
serious in the relevant region of small hole doping. We show that the difference in accuracy for 
different doping makes it very difficult to probe the phase separation instability using only energy 
calculations, even in the weak-coupling limit (U = At) where reliable results are available. By 
contrast, the knowledge of the charge correlation functions allows us to provide clear evidence of a 
spatially homogeneous ground state up to U = lOt. 



71.10.Fd, 71.45. Lr, 74.20.-z 
I. INTRODUCTION 

Since the discovery of high-temperature.-superconduc- 
tivity, the two-dimensional Hubbard modelEJ has been the 
subject of a huge amount of work. Indeed, it is widely 
accepted that, despite the complicate structure of these 
materials, a key role is played by the electron correla- 
tion in the CuO planes, well represented by the Hubbard 
model. 

This model and the related t — J model reproduce, 
at least qualitatively, some of the physical properties of 
cuprates. For example, at half- filling (i.e., when the 
number of electrons is equal to the number of sites) 
the ground state (GS) is an antiferromagnetic insulator, 
and upon doping the antiferromagnetism is strongly sup- 
pressed. The possible insurgence of superconductivity 
upon doping and the symmetry of the order parameter 
in the model are still open questions. Another important 
point to study is the stability with respect to phase sep- 
aration (PS) or charge-density waves (CDW) near half- 
filling. Indeed, both PS and CDW ha.ve been experi- 
mentally-observed in different systemsElB and different 
authorsoH have pointed out a possible relation between 
charge instability and superconductivity. 

Because of the strongly interacting nature of these sys- 
tems, important insights came from nonperturbative nu- 
merical methods. 

The exact diagonalization methods, e.g., Lanczos, are 
strongly limited by the exponential growth of the Hilbcrt 
space, and in practice it is possible to diagonalize only up 
to about 20 sites. A remarkable development of the ex- 
act diagonalization strategy is the density matrix renor- 
malization groupu, which allows us to compute the GS 
using an iteratively improved basis. Despite the accu- 
racy of this technique for one-dimensional and quasi-one- 
dimensional systems, there is no straightforward gener- 
alization for higher dimensions. 

Attractive alternatives are stochastic methods such as 
Quantum Monte Carlo (QMC) that do not suffer from 
severe limitations in terms of lattice size and hence al- 
low us to study fairly large systems. The simplest QMC 
simulation is the variational one, in which it is possible 



to compute the expectation values of the energy and all 
the correlation functions for a given wave function. The 
limitation of this method is that one has to guess a priori 
the form of the GS. 

QMC methods allow us to sample directly the GS by 
using projection techniques, the simplest one being the 
so-called power method, namely the simple idea that an 
iterative application of the Hamiltonian H filters out the 
ground state starting from an arbitrary state nonorthog- 
onal to it. These approaches suffer, in fermionic systems, 
from the so-called sign problem. Indeed, due to the an- 
tisymmetry of the wave function under permutations of 
two particles, one gets after a few power iterations oppo- 
site contributions leading to the cancellation of large and 
fluctuating weights. 

An alternative approach, which is quite efficient for 
the Hubbard model, is the auxiliary field quantum Monte 
Carlo (AFQMC)Q. The Coulomb interaction is linearized 
by using a Hubbard-Stratonovic transformation. In this 
way it is possible to apply exactly the operator e~@ H to 
a Slater determinant. In practice, AFQMC is very ac- 
curate only for small values of the Coulomb repulsion, 
whereas, when the interaction strength becomes compa- 
rable with the bare bandwidth, large fluctuations make 
the simulations rather unstable and inefficient, because 
the sign problem becomes particularly severe for strong 
Coulomb repulsion. 

A sophisticated implementation of the power method 
is the Green-function Monte Carlo (GFMC) applied to 
lattice models]. In principle, this method gives unbiased 
results but in practice, for fermion systems, the sign prob- 
lem makes any simulation prohibitive, because of statis- 
tical errors. Recently, many attempts have been made-io 
deal with fermions and overcome the sign problemBEj. 
The fixed-node approximation!!^ is introduced for the 
GFMC and replaces the original Hamiltonian with an 
effective one free of the sign problem, in which the nodes 
are fixed to be those of the so-called guiding wave func- 
tion. In practice, the most common choice is the best 
variational wave function available, i.e., the one with the 
lowest energy. The fixed-node energy is an upper bound 
to the true GS energy and the method is variational. 
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Moreover, it becomes exact if the guiding function is cho- 
sen to be the exact GS function. Although this approxi- 
mation can be uncontrolled and needs an external input, 
it is quite accurate with respect to both energy and cor- 
relation functions for small sizes. 

An alternative approach, recently proposed by one of 
us, is to use the GFMC with a stochastic reconfigura- 
tion (SR)B The main idea of this method is to use a 
suitable reference dynamics free of the sign problem to 
constrain the true dynamics. The fixed-node and a varia- 
tional dynamics are natural candidates as references. At 
each reconfiguration the amplitudes of the GS are chosen 
as small perturbations of the reference wave function am- 
plitudes such that the mixed averages of a few physically 
relevant quantities are conserved. The method is exact 
if all the operators in the Hilbert space are reconfigured, 
but it has been shown that a ipajor role is played by a 
very small number of operatorstn 

In this paper, we present a systematic study on the 
Hubbard model comparing different QMC techniques by 
studying in particular the stability of PS and CDW near 
half-filling for various values of the Coulomb repulsion. In 
Sec. H we introduce the Hubbard model and we discuss 



the implementation of our approximations, in Sec. [II we 
report our results, and finally in Sec. IV we give a brief 
summary. 



exp{— gj^i n i] n ii)i where g is a variational parameter, 
and J' is a Jastrow factor J — exp( ^ Y^i j v i,j^i^j )' 
where 7 is another variational parameter and Vij is taken 
from spin-waves theory£3. Care must be taken in the 
choice of the orbitals appearing in the Slater determi- 
nant. The most common choice is to take the orbitals 
from a Hartree-Fock (HF) approximation of the Hamil- 
tonian breaking the SU(2) spin rotation symmetry along 
the z axis. 
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where 



{m) = {mj) + {ni,i), 



(3) 



(4) 
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We consider only fillings which are closed shells for U = 
and where a solution with constant density (rij) = -j- = n 
and staggered magnetization (mj) = m are found. In this 
case, the HF many-body wave function can be written as 
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n <t n 
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II. THE MODEL 

We consider the Hubbard model on a square lattice of 
L sites with N = JVj + N± particles and Aj = Aj,, where 
Aj (Aj) is the number of spin-up (-down) particles. In 
order to study PS as close as possible to half-filling us- 
ing only closed-shell configurations, we consider square 
lattices tilted by 45° with L 2 = 2l 2 and I odd. In this 
way half-filling is a closed shell and the first doped closed 
shell has eight holes independent of L. The Hamiltonian 
reads 

H = -t ^ 4,<T C 3,<r + U ^2 n hl n i,l> (!) 

where ( ) stands for nearest neighbors, Cj )(T (c\ a ) destroys 
(creates) an electron with spin a at site i, and nj lCT = 
c \ a c i,a-- In the following, all energies are measured in 
units of t. 

As already pointed out above, in the presence of the 
sign problem the choice of the guiding wave function is 
crucial. Our wave function reads 



|*g) =Vs^oP g J\V), 



(2) 



where \T>) is a Slater determinant in which the orbitals 
are suitably chosen (see below), Vs z =o is the projector 
onto the subspace with A| = Aj, i.e., with zero total 
spin component along the z axis, V g is a Gutzwiller 
operator that inhibits the double occupancies V g = 



where the quasiparticles have definite momentum modulo 
Q = (tt, 7r) and definite spin, since the antifcrromagnetic 
order parameter is along the z axis 



Pi 



(7) 



k is in the reduced magnetic Brillouin zone, and Uk and 
Vk are defined in RefJlJ. It is worth noting that for U/t — > 
00, Uk = Vk = ;4|, namely the spin up and the spin down 

are in different sub-lattices (classical Neel state). 

In a previous workli3 the wave function (||) with \T>) 
given by Eq. (^]) has been found to be a rather poor 
approximation for large U/t at half- filling. In particular, 
in this representation the Jastrow factor J does not play 
any important role. 

We propose a wave function which is a straightforward 
generalization of the one successfully used for the Heisen- 
berg mode£3. The fundamental ingredient is to allow 
spin fluctuations perpendicular to the staggered magneti- 
zation. An easy implementation of this idea is to put the 
magnetization in the x — y plane allowing transverse fluc.- 
tuations along the z axis through a Jastrow-like factorEI 
This is achieved by a ^ rotation U y (w) around the y axis 
of the canonical operators: 



^(|)4^(f) = ^(4 T -4 



(8) 
(9) 
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The fermionic part of our guiding wave function is there- 
fore defined as a Slater determinant of the transformed 
orbitals, 



namely it is given by 



l,...,Nt 



1,...,N, 



\ v )= n 4,+ n ^,-i°>- 



(10) 

(ii) 

(12) 



Remarkably for U/t — > oo and half-filling, by construc- 
tion Eq. ( |l2|) becomes the Neel state with spin quantiza- 
tion parallel to the x axis, i.e. it has the correct Marshall 
sign on each of the 2 L configurations sampled by GFMC. 
In this limit, it is also clear why the Jastrow factor may 
be much more effective: being defined along the z axis, 
it allows us to sample the quantum fluctuation perpen- 
dicular to the staggered magnetization. In the previous 
case instead both the Jastrow quantization axis and the 
order parameter were parallel, and for U/t — > oo there 
is no way to sample any fluctuation, the only possible 
configuration being the classical one. 

It is worth noting that, although for the t — J model 
the d-wave BCS wave function with a spin-rotationally 
invariant density-density Jastrow factor represents a very 
accurate variational state, for the Hubbard model at 
small and intermediate coupling (U < lOt) the best 
choice for the variational and guiding wave function is 
given by a Jastrow-Slater determinant with rotated or- 
bitals Eqs. ([!(]) and (p"l"|). Indeed, although the d-wave 
BCS wave function is a singlet and does not break the 
SU(2) symmetry, it has a very poor variational energy for 
the Hubbard model. The quality of the variational energy 
obtained with our Jastrow-Slater determinant remains 
considerably better than the BCS one at half-filling and 
U < lOt even when the accuracy of the approximation 
is improved by the GFMC. Instead, in the doped case, 
the BCS wave function with GFMC is only slightly worse 
than the corresponding Jastrow-Slater determinant pro- 
posed in this work. This may suggest that antiferromag- 
netism is already suppressed at small finite doping, and 
d-wave superconductivity is a possible stable phase espe- 
cially at large U/t. 

An important systematic improvement of the wave 
function can be achieved by performing exactly one Lanc- 
zos step starting from |^g)j 



= (l+aJT)|* G >, 



(13) 



with a free parameter chosen to minimize the energy. 
This technique has been successfully used for the-£ — J 
model both to improve the variationaLcalculatiorJill and 
as a starting point for power methodalj. Henceforth, we 
will denote by VMC and LS the results obtained with 
the wave function Eqs. (pi) and (O), respectively, using 



variational Monte Carlo. Analogously, the symbols FN 
and FNLS will indicate the fixed-node approximation ap- 
plied to the wave function Eqs. (||) and ([l3]), respectively. 
Finally, the symbols SR will denote the stochastic recon- 
figuration approximation applied to the wave function 

In Tables H and O we report the energies of 18 and 
10 electrons on 18 sites, respectively. At half-filling we 
compare the results using Eq. (Q) with the Slater de- 
terminant | V) given by Eqs. (^) and ( |l2] ) for different 
approximations and values of U/t. Using Eq. (jl^) we 
obtain a sizeable improvement for large t/'s (U > lOt). 
Notice that for U = 20t the best variational result (the 
FNLS) with Eq. (|) is worse than the simple VMC with 
Eq. ("JJ). For 10 electrons the two Slater determinants 
give the same results. Indeed, for this doping the antifer- 
romagnetic order is strongly suppressed and the Jastrow 
factor J does not play any important role. 



III. RESULTS 

One of the most debated issues in strongly correlated 
electron models is the nature of thc-chArge distribution 
in their GS. Recently, many authorsc3~E3 have addressed 
the question of PS in the t — J model. It is well accepted 
that for J ^> t, holes tend to cluster together leaving the 
rest of the system in an antiferromagnetic state without 
holes. Although most of the calculations lead to the con- 
clusion that a critical value of J below which the GS is 
homogeneous exists, it is not clear what this value is at 
low doping, ranging from 0.5i and l.2t. Moreover, some 
authorsE3 have suggested that just before PS, the GS 
has charge modulations. Recently, an accurate numeri- 
cal study of a few chains on the Hubbard modelEZI has 
shown evidence of stripes, i.e., CDW oscillations. How- 
ever, this result appears limited to quasi-one-dimcnsional 
geometry, as also suggested by the authors. 

It is well known that PS is characterized by an infinite 
compressibility in the thermodynamic limit. The com- 
pressibility can be related to the curvature of the energy 
with respect to the electron density, 



X 



d 2 E 



\ dn 2 



(14) 



From the above definition we have that a divergent x 
corresponds to a vanishing curvature of the energy as a 
function of density. Therefore, PS can in principle be 
detected by means of energy measurements for various 
densities. This is an appealing property, since the accu- 
racy on the energy is usually better than that of any other 
observable for most numerical methods. Many previous 
numerical studies of PS have therefore concentrated on 
the energy curve, or equivalently, on the energy per hole 
eh{S) = [e(8) — eH]/S, where e(S) is the energy per site at 
a hole density S = 1 — n and en — e(0) is the energy at 
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half-fillingpl. If ^ diverges, e^iS) is flat in the thermody- 
namic limit and develops a minimum for 8 = S c in finite 
systems, due to the finite positive surface energy at the 
phase boundary. .— . 

This approach has been pursued by Cosentini et al&B, 
using FN calculations. They found that there is a large 
region of PS in the phase diagram, at least for U > lOt, 
clearly in contrast with what is found in the t — J model. 
In this paper, we consider the Hubbard model and we 
show that a study of PS instability is very difficult using 
only energy calculation. Instead, a careful calculation of 
charge-correlation functions strongly indicates that the 
GS is homogeneous. 

In order to show that the energy calculations may over- 
estimate the tendency to a PS instability, it is important 
to compare the GFMC results with some exact reference 
results. Previous studies have shown that it is important 
to consider relatively large lattice sizes since finite-size 
effects favor PSt3. We need, therefore, a reference re- 
sult for large lattices, where exact diagonalization is not 
available. In the case of the Hubbard model for small U, 
the AFQMC is almost exact and represents the reference 
we need. As stated in the preceding section, we consider 
only closed shell doping because, at least for small U, 
huge finite size effects affect the physical properties in a 
rather drastic way. For instance, the large bare density 
of states near half-filling determines an. unphysical and 
spurious PS up to the first closed shellc3. 

In Fig. [I], we show the accuracy of the GFMC results 
obtained with different approximations compared with 
the AFQMC ones for a 162-site lattice and U = At. For 
this coupling value, AFQMC does not provide evidence 
for PS. We plot \E{§) - E hs (S)]/E hs (8), where E(S) and 
E hs (6) are the energies of GFMC and AFQMC, respec- 
tively, for a doping 5. Besides the improvement in the 
absolute accuracy, the curves get flatter and flatter im- 
proving the approximation, but only the SR accuracy is 
almost doping- independent. In other words, we need a 
very accurate calculation to eliminate the spurious de- 
pendence of the variational energy upon dopingc3. Even 
for the best variational method, the FNLS, although the 
accuracy on the energy is for all doping less than 1%, 
the difference in accuracy between, for example, the half- 
filled case and the first closed shell is still sizable. This 
difference is very important, because it represents just 
the energy scale determining or ruling out PS. 

In Fig. H the function eh{S) is shown for the FN, FNLS, 
SR, and AFQMC methods. We need to use SR to exclude 
the occurrence of PS, where even the FNLS data would 
imply PS. The reason for this disappointing situation is 
that all the known variational methods are still too de- 
pendent on the guiding wave function. With the previous 
analysis, the resolution in energy necessary to detect or 
rule out PS is very hard to reach with statistical methods, 
especially for large U/t. 

On the other hand, GFMC methods have proven to 
be reliable not only for energy calculations but also for 
correlation functions such as N(q) — (n q n^ q ), where n q 



is the Fourier transform of the electron densityc-l For a 
phase-separated system, there are strong fluctuations in 
the density for small momenta and N(q — > 0) is expected 
to be strongly enhanced for small momenta, that is, for 
\q\ ~ ^r, where £ is the characteristic length of the phase- 
separated region. Moreover, if \ diverges, also N(q — > 0) 
diverges, yielding an alternative tool to probe PS. 

This method turns out to be more reliable, since it is 
based on a single calculation for a given doping value, 
whereas the evaluation using eh(S) involves a compari- 
son between energies obtained by different simulations 
for different fillings, with the corresponding guiding wave 
function having different accuracies. Indeed with GFMC, 
N(q) has been proved to be a very sensitive tool to look 
at for detecting PS. In the t — J model, N(q) has a 
very different shape for stable and unstable systems. A 
clear peak at the smallest q indicates PS even when, as 
shown in tho— inset of Fig. ||, J/t is very close to the 
PS boundarye.3. Moreover, from N(q) it is also possi- 
ble to extract information about .charge fluctuations at 
finite g's, related to CDW. In RefEJ we have shown that 
in the t — J model for J — OAt, N(q) has some peaks 
at finite g's, surnxisingly near to what was found in re- 
cent experimentsEII. The knowledge of N(q) allows us to 
extract more general results with respect to the simple 
study of eh(S). Furthermore, N(q) is found to be much 
less size dependent than eh(S). For the t — J model there 
is no appreciable difference between 98-site and 162-site 
lattices and for J = OAt. Both calculations suggest that 
there is no PS, whereas an analysis using the FN approx- 
imation of e?j, (6) should lead to PS for 98 sites and to a 
homogeneous state for 162 sites. Indeed, in this case PS 
is only a size effect and a homogeneous state is found by 
increasing the accuracy of the method or by increasing 
the lattice size. 

We computed N(q) by means of the forward- walking 
techniquetj, within the FNLS approximation, at half- 
filling and for the first few closed-shell configurations on a 
162- and a 98-site lattice. The evaluation of the density- 
density correlation function is in principle possible even 
within SR by numerical differentiation of the energy with 
respect to an external field coupled to N(q). However, 
this approach is very demanding and does not give a 
significant improvement on the FNLS results, which are 
very accurate. Indeed for the smallest q vector for 90 
electrons on the 98-site lattice, i.e., q = (27r/7, 2n/7), we 
found the the AFQMC gives N(q) = 0.0932(2) and the 
FNLS gives N(q) = 0.096(1). 

In Fig. ||, N(q) is shown for U = At at half- filling for 
a 162-site and a 98-site lattice and for 154 electrons on a 
162-site lattice and 90 electrons on a 98-site lattice. No 
sign of divergence, and consequently of PS or CDW, is 
seen in the data. We also notice that the two sets of 
points for the half-filled systems lie on the same curve, 
showing that we have substantially reached the thermo- 
dynamic limit. For this value of U, N(q) is essentially 
featureless for all doping we considered, suggesting that 
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there are no charge instabilities at any finite length. The 
smallest doping we considered is S ~ 0.049 and we cannot 
exclude that for smaller doping PS or CDW are present. 

In order to investigate smaller doping, we should con- 
sider larger lattices. Unfortunately, the accuracy of the 
approximations considered decreases when increasing the 
size of the system and the 162-site lattice represents the 
largest lattice where the accuracy is acceptable. In Fig. ||, 
we report [E(0) — Eh s (0)]/Eh s (0) for various sizes and 
for different approximations: from the 18 sites to the 
162 sites, the accuracy of FNLS changes from less than 
0.1% to about 0.5%. These indications prevent us from 
considering sizes larger than the ones presented in this 
paper. 

Now we turn to larger Coulomb interactions and con- 
sider U = lOt, where the AFQMC results are not re- 
liable due to large fluctuations. In principle, GFMC 
techniques do not suffer from intrinsic limitations in the 
large-coupling regime and it is possible to consider any 
value of U. In practice we need an accurate knowledge of 
the nodes, i.e., an accurate guiding wave function. Our 
choice, Eq. (||), with orbitals given by Eqs. (|H]) and @ 
is a very good approximation for the half-filling case. In 
Table [j], we report the energies for various methods for 
18 electrons on 18 sites at different U's. Although all 
the approximations are quite size-dependent, the wave 
function becomes more and more accurate by increasing 
the Coulomb potential. Therefore, we expect that it also 
gives a good starting point at least close to half-filling. 

We present results for U — lOt, for which previous 
FN calculations based on eh (6) and a less accurate wave 
functionE£l have shown PS up to 5 ~ 0.15. Indeed, if 
we use eh(5) as a probe for PS, we find that the phase 
diagram shows-a large instability region, confirming the 
results of Ref.EL As for the U = At case, this instability 
is very likely to be a spurious effect, a consequence of 
the different energy accuracy for different doping. This 
possibility, which cannot be proved without knowing the 
exact energies at strong coupling (at present impossible), 
is instead very clearly supported by the calculation of the 
charge-correlation functions. 

Figure || displays N(q) for the same fillings of Fig. || 
and for 138 electrons on the 162-site lattice, which cor- 
responds to 8 ~ 0.148. All the correlation functions are 
definitely nondivergent for q — > and are qualitatively 
similar to the U = At case, indicating that the system is 
far away from a PS instability. Furthermore N(q) does 
not show peaks at any finite momenta for this Coulomb 
interaction. This finding shows that the Hubbard and the 
t — J model may have different behaviors as far as charge 
correlations are concerned. Indeed, by diagonalizing ex- 
actly the 18-site lattice, we find that charge fluctuations 
for the Hubbard and for the t — J model are quite differ- 
ent in the spall doping region for U = lOf and J = 0.4f, 
respectively^! 



IV. CONCLUSIONS 

An extensive GFMC analysis of the Hubbard model 
at low hole doping has been carried out. In particu- 
lar, we have focused on the possible instability of the 
model with respect to PS and CDW. Comparing GFMC 
results with AFQMC ones in the weak-coupling region 
(U = At), we show that detecting PS by means of energy 
results requires a very accurate calculation at all elec- 
tronic densities. Indeed, the accuracy of the energy is 
strongly dependent on the electron density, and the sig- 
nature of PS based only on energy calculations is clearly 
affected by this bias, leading to a spurious region of PS 
instability. In the case of the Hubbard model, this is 
particularly relevant because, while we are able to give a 
very good description of the half-filled case, in which the 
GS is an antiferromagnetic insulator, we are not aware of 
equally accurate descriptions of the doped state. Even for 
U = At, it is necessary to use the really accurate SR and 
AFQMC technique to eliminate the doping dependence 
of the accuracy and to rule out PS. 

On the other hand, PS (and CDW) instability can be 
probed more easily using charge correlation functions. 
This approach has various advantages. First, it is found 
that N(q) has very small size effects and the thermody- 
namic limit is reached with about 100 sites, both for the 
Hubbard and the t — J model. Second, the information 
contained in N(q) does not depend on different densities, 
implying that a different accuracy as a function of doping 
does not introduce any external bias. 

Instead of using energy calculations, which are very 
expensive at moderate and large U's, we calculate the 
charge correlation functions and we are able to find clear 
evidence for the absence of PS up to U = lOf in the low 
doping regime. 
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FIG. 1. Relative accuracy of various GFMC techniques 
with respect to AFQMC for a 162-site lattice with U = At 
as a function of filling 8. From top to bottom VMC (empty 
triangles), FN (empty hexagons), LS (empty squares), FNLS 
(empty circles), and SR (full squares). Lines are guides to the 
eye. 
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FIG. 3. N(q) for U = At, 162 electrons on 162 sites (empty 
circles), 98 electron on 98 sites (empty squares), 154 electrons 
on 162 sites (full circles) and 90 electrons on 98 sites (full 
squares) . Lines are guides to the eye and error bars are smaller 
than points. T = (0, 0), X = (tr, n), M = (ir, 0). In the inset: 
N(q) for the t — J model, J = 0.6i, 156 electrons on 162 sites 
(squares), and 94 electrons on 98 sites (circles). 
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FIG. 2. Energy per hole eh(5) for a 162-site lattice with 
U = At. From top to bottom FN (empty hexagons), FNLS 
(empty circles), SR (full squares) , and AFQMC (full circles). 
Lines are guides to the eye. 
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FIG. 4. Relative accuracy of various GFMC techniques 
with respect to AFQMC for different lattices (L = 18, 98, 162) 
and U = At. From top to bottom: VMC (empty triangles), 
FN (empty hexagons), LS (empty squares), and FNLS (empty 
circles). Lines are guides to the eye. 
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FIG. 5. N(q) for E7 = lOt, 162 electrons on 162 sites 
(empty circles), 98 electrons on 98 sites (empty squares), 154 
electrons on 162 sites (full circles), 90 electrons on 98 sites 
(full squares), and 138 electrons on 162 sites (full triangles). 
Lines are guides to the eye and error bars are smaller than 
points. T = (0,0), X = (ir,n), M = (tt,0). 
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-0.4034(1) 
-0.2023(1) 
-0.9460(1) 
-0.4382(1) 
-0.2293(1) 



-0.9514(1) 
-0.4284(1) 
-0.2195(1) 
-0.9547(1) 
-0.4451(1) 
-0.2232(1) 



-0.9520(1) 
-0.4154(1) 
-0.2060(1) 
-0.9553(1) 
-0.4428(1) 
-0.2310(1) 



-0.9556(1) 
-0.4316(1) 
-0.2225(1) 
-0.9576(1) 
-0.4470(1) 
-0.2337(1) 



TABLE I. GS energies for 18 electrons on 18 sites as a 
function of U/t using H and (H) as Slater determinant. 



U E e 



Evmc 



Efn 



E l < 



Efnls 



4t -1.1299 -1.1124(1) -1.1218(1) -1.1229(1) -1.1263(1) 
lOt -1.0193 -0.9749(1) -1.0006(1) -0.9997(1) -1.0098(1) 
20t -0.9598 -0.8983(1) -0.9354(1) -0.9253(1) -0.9450(1) 

TABLE II. GS energies for 10 electrons on 18 sites as a 
function of U/t. 
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